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Abstract 

Two conditions on symmetries are identified as necessary for a general linear scattering system to 
be able to rotate the linear polarisation of light: Lack of at least one mirror plane of symmetry 
and electromagnetic duality symmetry, which is equivalent to the conservation of the helicity of 
light. When the system is a solution of a single species of particles, the lack of at least one mirror 
symmetry leads to the familiar requirement of chirality of the individual particle, and, according to 
the analytical and numerical evidence presented in this paper, the solution preserves helicity if and 
only if the individual particle itself preserves helicity. However, in forward scattering the helicity 
preservation condition on the particle is relaxed: We show that the effective rotational symmetry 
of such system leads to helicity preservation in the forward scattering direction independently of 
any property of the particle. These results advance the current understanding of the phenomena 
of molecular optical activity and provide insight for the design of polarisation control devices at 
the nanoscale. 



1 



An object which cannot be superimposed onto its mirror image is said to be chiral. Chi- 
rality is entrenched in nature. For instance, some interactions among fundamental particles 
are not equivalent to their mirrored versions Jj. Also, the DNA molecule is a double helix 
strand of chiral amino acids. The understanding and control of chirality has become impor- 
tant in many scientific disciplines. In chemistry, the control of the chiral phase (left or right) 
of the end product of a reaction is crucial, since the two versions can have very different 
properties. In nanoscience and nanotechnology, chirality plays an increasingly important 
role 

Electromagnetic waves can also be chiral. The photon has its chiral character mapped 
onto a binary property, the helicity {4], From the considerations above, and since 
electromagnetic waves are used to interact with matter at the nano, molecular and atomic 
scales, it is not surprising that the interaction between chiral light and chiral matter is an 
important subject of study. Interestingly, the subject is quite old and, from the beginning, 
has always been associated with the rotation of the linear polarization of light. For instance, 
Biot discovered that when lightpropagates through a solution of certain types of molecules, 

6] . Commonly referred to as molecular optical activity, the 



its linear polarization rotates 

study of its root causes has a long history reviewed for instance in [?], 



and 



9|. In 1848, 



asteur identified the lack of any mirror symmetry of the molecule as a necessary condition 



10| . He called it "dissymetrie moleculaire" and by it Pasteur meant non-superimposability of 



the molecule and its mirror image: Chirality. Nowadays, this necessary condition is assumed 



to also be sufficient, and the exceptions to the rule are explained by other means [8], 11]. 
A comprehensive theoretical study of optical activity based in symmetry principles can be 
found in |9(, and the modern theoretical and computational methods for optical activity 
calculations are reviewed in 12| . Current investigations of optical activity in metamaterials 
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15] are aimed at obtaining polarization manipulation devices suitable for integrated 



nanophotonics. 

In this article we rigorously study the symmetries of a general linear system that rotates 
linear polarization. We find the lack of at least one mirror symmetry and invariance under 
generalized electromagnetic duality transformations to be necessary symmetry conditions for 
such system, generalized duality invariance is equivalent to the preservation of the helicity 



of light 
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17| . We then apply these general results to a mixture of randomly oriented 



replicas of a single particle with the aim of identifying the restrictions that each of the two 



necessary conditions imposes on the electromagnetic response of the individual particle. We 
find that the lack of at least one mirror symmetry of the mixture translates into the condition 
of "dissymetrie" for the individual particle. With respect to duality symmetry, we provide 
analytical and numerical evidence which strongly suggests that, in general, the necessary 
and sufficient condition for the mixture to preserve helicity is that the individual particle 
preserves helicity. Our results can be reconciled with the current understanding of molecular 
optical activity by restricting the measurements to the forward scattering direction only. 
Due to the large number of randomly oriented particles, the mixture acquires an effective 
rotational symmetry, which is shown to lead to the conservation of helicity in the forward 
scattering direction independently of any property of the individual particle. Therefore, in 
the forward scattering direction, a solution of a chiral molecule can exhibit optical activity 
without the individual molecule having to preserve helicity. In an ordered system or in any 
non-forward scattering direction, explicit helicity preservation by the individual particle is 
needed for optical activity. 




FIG. 1. (Color online) Transformation of input linear polarization states (diagrams on the left) to 
rotated elliptical output polarization states (diagrams on the right). The rotation angle (3 and the 
major to minor axis ratio (ellipticity) of the ellipse do not depend on the input angle a. 

Throughout the paper, we will assume monochromatic electromagnetic fields with a har- 
monic exp(— iut) time dependence. Consider the polarization transformation illustrated by 
FigfTJ A linearly polarized input transforms into an elliptically polarized output with its 
major ellipse axis rotated by a fixed angular quantity (3 with respect to the angle a of the 
input linear polarization. Both the rotation angle /3 and the major to minor axis ratio 
(ellipticity) of the ellipse are independent of a. This is a straightforward generalization 
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of the experimentally observed molecular optical rotations. However, other definitions are 
sometimes used in the metamaterial literature where f3 depends on a jlol \l8\. Our aim is 
to identify what restrictions this class of transformations imposes over a general conversion. 
This is most easily achieved expressing the input linear polarization state [cos (a) , sin(a)] T 
in the circular polarization basis 
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where a general transformation of the input reads 
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The angle of the major ellipse axis with respect to the horizontal axis is 6 — | arg (EiE r } . 
According to our specification 29 = 2(a + /3) V a, which then forces EiE r = r]exp(i2(a + /?)) 
where 77 is a real number. Using ([1]) we then obtain the relationship: 



EiE T 



ac 



ad* exp(i2a) + be* exp(—i2a) + bd* — rjexp(i2(a + (3)), 



which is valid for all a and hence imposes b = c = and gives 2(3 = 2(arga — argo?). 
The most general matrix which meets the requirements is hence diagonal. We conclude 
that our specified transformation is equivalent to the conservation of circular polarization 
states. From now on, we will refer to such a transformation as a generalized rotation of 
linear polarization, and abbreviate it by GRLP. 

We now consider an electromagnetic scattering situation, where an incident field interacts 
with a general linear scatterer. As a result of the interaction a scattered field is produced. 
Both incident and scattered fields can be decomposed in plane waves. We impose that the 
relationship between the polarizations of any pair of incident and scattered plane waves is of 
the kind depicted in FigJT] which results in a transformation matrix of the required diagonal 
form. The conclusion is that the system preserves the circular polarization of any plane wave, 
which is equivalent to the preservation of the electromagnetic helicity. Therefore, the system 
is invariant under the transformations generated by helicity: generalized electromagnetic 
duality HQ- 

In the case of the photon, helicity can take the values ±1. It is defined as the operator 
resulting from the projection of the angular momentum onto the direction of the linear 



momentum J4], i.e. A = J • P/|P|. For a given plane wave of momentum vector p, the 

two possible states of definite angular momentum along p (J p ) correspond to the states of 

circular polarization of the plane wave with eigenvalues equal to ±1. 

□ 

In the same way that angular momentum generates rotations [19j, helicity is the generator 
of the generalized duality transformations [20| : 



E' = cos 7E — sin 7H, H' = sin 7E + cos 7H, 



(2) 



where 7 is a real angle. The typical duality transformation E — )• H, H — > — E is recovered 
for 7 = —tt/2. Helicity preservation and invariance under the transformation in fl2]) are 
hence equivalent conditions. The GRLP condition imposed on our system has lead us 
to conclude that its electromagnetic response must be invariant under generalized duality 
transformations. 

To proceed with the study of the symmetries of our system, we introduce a concise 
notation that expresses the action of the system 5* as a linear operator which takes input 
plane waves of a given momentum p and helicity ±, |p, ±) into output plane waves |p, ±). 



S = J dp J dp (ag[p, +) <+, p| + dg|p, -) <-, p|) . 



(3) 



The orthogonality relationships (A,p|p', A') = 0a-a"5(p — p') are assumed. The absence 
of cross-terms |p, — )(+, p|, |p, +)(— , p| in <^ reflects the helicity preservation condition. 
Taking a pair (p, p), the 2x2 matrix that specifies the conversion is 
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|dP| exp(i arg rfg) 



The dependence of the ag and dg on the momenta (p, p) allows for different transformations 
of the required diagonal kind with rotation angles = |(argaP — argd^). 

We now ask the following question. Let us assume that there exists a pair (p, p) for which 
/?p 7^ 0. What can be said about the symmetries of the system? To answer this question, 
we consider the mirror operation M pp across the plane defined by (p, p) and assume that 
the system possesses this mirror symmetry: M~pSM p p = S. This mirror operation leaves 
the momentum vectors invariant and since any mirror operation flips the helicity value, the 
states transform as |p, =f) = M P p|p, ±), |p, =1=) = M p p|p, ±). Using them, and that the 
mirror operator is unitary: M~p = M pp we can see that, if the system is invariant under 



this mirror transformation, 



4 = (+,p|£1p,+) = (+,p\Ml P SM pP \p,+) 

= (-,p|5|p,-) = 4^/3P = 0. 

Therefore, if there exists a pair (p, p) for which /3£ 7^ 0, what we can say about the system 
is that it does not have that particular mirror symmetry: M~~SM p p 7^ S. Summarizing the 
discussion up to this point: By imposing that the system performs GRLP for any pair of 
input output plane waves we have found that it must preserve helicity and if at least one of 
the rotation angles is to be different from zero, it must lack at least one mirror symmetry. 
These are two necessary conditions for GRLP. Expressed with the help of commutators 21] 
between operators and denoting as the mirror operation across the plane perpendicular 
to vector a, our findings are: 

GRLP [S, A] = and 3 a s.t. [3, M & ] ^ 0. (4) 

It can be shown that the two necessary conditions are not sufficient. It can also be shown 
that, when further assuming rotational symmetry, the three conditions together are indeed 
sufficient for nonzero GRLP. 

Having derived for a general linear system, we now turn to the study of GRLP by a 
mixture containing a large number of randomly oriented scattering particles immersed in an 
isotropic and homogeneous medium. We assume that the mixture has a linear response and 
contains only one kind of particle. Our aim is to investigate what conditions (jl]) imposes on 
the scattering operator S u of the individual particle. For this, we make use of the theory 
of independent random scattering to compute the Mueller matrix of the mixture. The 
Mueller matrix relates the input Stokes parameters with the output Stokes parameters. The 



theory of independent random scattering 22| 



231 ] tells us that, if the individual particles 



are sufficiently separated and the number of particles tends to infinity, the Mueller matrix 
of the total solution Lg(p, p) can be computed as the average sum of the Mueller matrices 
for all possible orientations of the individual particle. If /(•) is the function that converts a 



24j, we have that 



2x2 scattering matrix to its corresponding Mueller matrix 

L s (p,p)=n j dRf(S*(p,p)). (5) 

where n is the density of particles per unit volume, j dR indicates the sum over all possible 
rotations and 5^(p, p) is the 2x2 scattering matrix with coefficients (A, pl-R^S^-Rlp, A). It 



is important to note that due to the integral over all rotations, equation <^ is only exact 
in the limit of infinite number of particles. From now on, we will take (jSJ) as an effective 
response for the mixture and identify the results that explicitly rely on the J dR average. 

We start with the condition concerning mirror symmetry: 3 a such that [S, M&] 7^ 0. 
The Mueller matrix of the mirror system can be shown to be 

L Mt S M,(P>P) = n ° / dRf((\,p\rfMis u M & R\p,\)). (6) 

Lack of the mirror symmetry M& implies that Ls(p, p) 7^ ^mIsm- (p> p) ^ or a ^ l eas t one 
pair (p,p). 

Now, let us assume that the individual particle possesses a symmetry of the rotation 
reflection kind |3Jj: MsRs (g), with m a positive integer. When we assume any of these 
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symmetries for S u , the argument of /(■) in (jSJ) can be written as (A, p\R^ R} S u RR\p, A) 
Then: 

L Mt S M a (P'P) = n o / dRf({\,p\rfi&SjLR\p,\))=Ls(p,p). 

The second equality follows from the fact that, when R covers all possible rotations once, 
RR also covers all possible rotations once. This is an application of the re-arrangement 
lemma from group theory jj]. 

We have just proved that for the mixture to lack one mirror symmetry the individual 
particle must not have any symmetry of the type M&R& (f^j. The lack of all of these 



symmetries (for all m) is equivalent to the particle being chiral 11]. Some qualitative 
consideration suffices to realize that such condition on the particle is not only necessary, 
but also sufficient for the mixture to lack not one but all mirror symmetries. It is actually 
impossible for the random mixture to possess one mirror symmetry without possessing them 
all. 

We now turn to the duality condition [S, A] = 0. If we impose [5 U ,A] = 0, it is easy 
to show that L s (p, p) = L A ^ SA (p,p) for all (p,p). Importantly, the f dR averaging does 
not need to be invoked in such proof. This result also applies to an ordered mixture or 
a mixture with a small number of particles. We now ask: Could it be that due to the 
averaging and/or randomness of the mixture, [S u , A] 7^ but [S, A] = 0? Using again 
the re-arrangement lemma (which makes use of the j dR averaging) allows to show that 

(p? p) — L/r\sr{P^ p)j f° r an y rotation R. Such effective rotational symmetry ensures 
helicity preservation in the forward scattering direction (p = p): Helicity being the angular 



momentum along the momentum axis of the input and output momenta, it must be preserved 
by a rotationally symmetric system when p = p. Having used the average over all possible 
rotations f dR to derive this result, it will not apply to ordered systems or systems with a 
small number of particles. 

We now provide analytical and numerical evidence that helicity is not preserved when 
p ^ p unless [S u , A] = 0. The Mueller matrix L of a helicity preserving system should 



preserve the two input Stokes vectors of well defined defined helicity 
imposes relations between the matrix coefficients L^: 



1 ±1 



Ln — L44, L14 



L 



31 



L 



24 



L 



34 



0. 



which 



(7) 



In general, equation (JSJ) must be evaluated numerically. We later provide numerically 
obtained values for mixtures of spherical, conical and helical particles. For the simple case 
of small (w.r.t. the wavelenght) spherical particles with electric constant e and magnetic 
constant \i = 1, there is an analytical expression for S^p, p) [23| which allows to compute 
Ls(p,p) , 

6 - 1 



2„3 



47r(e + 2) 
cos 2 (x P p) + 1 cos 2 (x P p) - 1 
cos 2 (x P p)-l cos 2 (x P p) + l 





(8) 





2cos(x P p) 

2cos(x P p) 

where k is the wavenumber, a the radius of the sphere and Xpp is the angle between the input 
and output momentum vectors. For general x p p, matrix ([S]) violates the helicity preserving 
conditions ((71). We conclude that a necessary condition for a solution of small spheres to 
preserve helicity is that the individual sphere preserves helicity. We have already discussed 
that such condition is also sufficient. 

To investigate whether the conclusions reached for small spheres also hold for mixtures 
of other kinds of particles and sizes, we numerically computed the rotational average §5§ 
for small conical, and helical particles and for spheres of different sizes, with e = 2.25 and 
[i = 1 immersed in vacuum. From a recently obtained result 2(| regarding the conditions for 
duality symmetry (helicity preservation) on the macroscopic Maxwell's equations, we know 
that a particle of arbitrary shape with electric constant e and magnetic constant [i immersed 



S 



in a solvent (e s ,fi s ) would preserve helicity if and only if e//i = e s /// s . Then, particles with 
e = 2.25, fj, = 1 are non-dual when immersed in vacuum and hence do not preserve the 
helicity of light. 

To measure the degree of helicity transformation in each case we use the following metric 

on the resulting Mueller matrices: 

p_ (Ln + Lu ~ (L41 + L44)) 2 

2 (L n + Lu - (L41 + L44)) 2 + (L n + L u + (L41 + L44)? ^ 

(Ln — L14 + (L44 — L44)) 2 

2 (£u — L 14 + (L41 — L44)) 2 + (Ln — L14 — (L41 — L44)) 2 

The first (second) line in (Q is the relative helicity change effected by the Mueller matrix 

on a Stokes vector of well defined positive (negative) helicity. T = for a helicity preserving 

Mueller matrix ((7j), and T = 1 for a helicity flipping Mueller matrix. 



See [26j] for the details of the computations. Figj2]^a) plots r(x P p) for cones and helices 
of dimensions ~ A /6, where A is the wavelength. Figj2](b) plots T(x P p) for spheres of 
diameters ~ [A /2, A , 2A ]. All cases show that the solutions do not preserve helicity, 
strongly suggesting that, in general, for the solution to preserve helicity, the individual 
particles must preserve helicity. In all cases T(0) = in full agreement with our previous 
discussion about forward scattering (p = p). In all cases T(7r) = 1, indicating that for 
backward scattering (p = — p), helicity is always exactly flipped. This is indeed a general 
result due to the effective rotational symmetry of the solutions. The angular momentum 
along — p must be the same as the input angular momentum along p since they share the 
same axis. Since helicity is A = J-P/|P|, it must exactly change sign due to the preservation 
of J and the sign change in P in the backward scattering direction p = —p. 

The current understanding of the phenomena of molecular optical activity is that chi- 
rality of the individual molecule is the only necessary and sufficient condition. We can 
only reconcile our results with the current views by restricting the measurements to the 
forward scattering direction only. As we have discussed, the large number of randomly ori- 
ented particles endows the mixture with an effective rotational symmetry, which is shown 
to lead to the conservation of helicity in the forward scattering direction independently of 
any property of the individual particle. Therefore, in the forward scattering direction, a 
solution of a chiral molecule can exhibit optical activity without the individual molecule 
having to preserve helicity. In an ordered system or in any non-forward scattering direction 
the effective rotational symmetry disappears and explicit helicity preservation (generalized 




FIG. 2. (Color online) Helicity transformation metric T(xpp) results obtained from the numerical 
computation of the Mueller matrices of vacuum solutions of different kinds of particles with e = 2.25 
and fj, = 1. r = corresponds to helicity preservation and T = 1 to complete helicity change, 
r = for all Xpp would correspond to a helicity preserving solution, (a) Spheres, cones and helices 
of dimension ~ Ao/6. The wavelength of light was Ao = 632 nm. (b) Spheres with diameters 
~ [Ao/2, Ao, 2Ao]. All cases break helicity preservation. In the small particle case (a), the results 
are very similar (small differences not visible in the figure) independently of the kind of particle. 
This can be explained because all small particles can be treated in the dipolar approximation. 
When the size of the particles is comparable to the wavelength (b), we observe a more complex 
behavior of T(xpp) with oscillations with T close to zero for some angles. We speculate that this 
behavior is related to the excitation of higher (than dipole) multipolar moments. 

electromagnetic duality symmetry) by (of) the individual particle in the solvent is needed 
for optical activity. 

In this article, we have used the formalism of symmetries and conserved quantities to 
study a class of electromagnetic transformations which we have named generalized rotation 
of linear polarization (GRLP). We have identified two symmetry conditions necessary for 
a general electromagnetic linear scattering system to exhibit GRLP: Lack of at least one 
mirror symmetry and generalized duality symmetry (helicity preservation). For a random 
mixture of a single species of particle immersed in an isotropic and homogeneous medium, we 
have investigated the restrictions that the two necessary conditions impose on the individual 
scattering properties of the particle. We have proved that the individual particle must be 
chiral and we have provided strong analytical and numerical evidence indicating that the 
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individual particle in the solvent must itself have an electromagnetic response invariant 
under generalized duality transformations, that is, it must preserve the helicity of light. Our 
results advance the current understanding of the phenomenon of molecular optical activity. 



20] used in this article, 



The general conditions in (J3J) together with the main result of 
also provide insight that may assist the design of polarization control devices, particularly 
at the nanoscale where metamaterials are used to engineer effective electric and magnetic 
constants. 
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